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Abstract 

We study the optical absorption spectra in a polar ferrimagnet GaFeOs. We consider the El, 
El and Ml processes on Fe atoms. It is shown that the magnetoelectric effect on the absorption 
spectra arises from the El-Ml interference process through the hybridization between the 4p and 3d 
states in the noncentrosymmetry environment of Fe atoms. We perform a microscopic calculation 
of the spectra on a cluster model of Fe06 consisting of an octahedron of O atoms and an Fe atom 
displaced from the center with reasonable values for Coulomb interaction and hybridization. We 
obtain the magnetoelectric spectra, which depend on the direction of magnetization, as a function 
of photon energy in the optical region 1.0 — 2.5 eV, in agreement with the experiment. 

PACS numbers: 78.20.Ls, 78.20.Bh, 78.40.-q 
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I. INTRODUCTION 



It is known that the breaking of time-reversal symmetry in magnetic materials gives rise 
to interesting magneto-optical effects such as the double circular reflection for circularly po- 
larized light and the Faraday effect for linearly polarized light.- When the spatial inversion 
symmetry is further broken, for example, in polar or chiral materials, novel magneto-optical 
effects were expected to come out.- Those effects are known as the nonreciprocal directional 
dichroism or magnetochiral dichroism, and have been extensively studied .-^S^ Among a 
variety of compounds, Cr 2 03 is one of the most investigated compounds. The magnetoelec- 
tric effect, that is, a linear relation between the magnetic and electric fields in matter was 
proved in 1950s.->2 Later, the nonreciprocal rotation and ellipticity of light were measured,— 
and were successfully analyzed by using a ligand field model for Cr atoms.— 

Another notable compound is GaFeC>3, which was first synthesized by Remeika.— This 
compound exhibits simultaneously spontaneous electric polarization and magnetization at 
low temperatures. The large magnetoelectric effect was observed by Rado.— Recently, un- 
twinned large single crystals have been prepared,— and the optical absorption measurement 
has been carried out with changing the direction of magnetization.— It has been found that 
the absorption intensity in the region of photon energy 1.0 — 2.5 eV changes with reversing 
the direction of the magnetization. The purpose of this paper is to analyze in detail this 
phenomenon by carrying out a microscopic calculation of the spectra and to elucidate the 
microscopic origin. Although several qualitative arguments have been done,— ^ as far as we 
know, the spectra have not been calculated yet as a function of photon energy. 

The crystal of GaFeC>3 has an orthorhombic unit cell with the space group Pc2in.— 
The magnetic moments at Fel and Fe2 sites align antiferromagnetically along the ±c axis. 
The actual compound, however, behaves as a ferrimagnet,— which reason is inferred that 
the Fe occupation at Fel and Fe2 sites are slightly different from each other.— Each Fe 
atom is octahedrally surrounded by O atoms, and slightly displaced from the center of the 
octahedron; the shift is 0.26 A at Fel sites and — O.llA at Fe2 sites along the b axis.— 
Thereby the spontaneous electric polarization is generated along the b axis. We neglect 
slight distortion of octahedrons, since their contributions are expected to be small to the 
El-Ml terms. There are two kinds of octahedrons with respect to the direction of Fe shift, 
as illustrated in Fig. [D 
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FIG. 1: Two kinds of octahedrons of oxygen atoms (white circles). Fe atoms (black circles) are 
displaced from the center of the octahedron O to the off-center O' along the b axis by amount 8; 
5 = 0.26A at Fel sites and 5 = -O.llA at Fe2 sites. 

In the analysis of optical absorption, we assume that the photon propagates along the a 
axis in accordance with the experimental situation.— Restricting the processes only on Fe 
atoms, we derive the explicit forms of El, E2, and Ml transitions. We find that the E2 
transition matrix elements are much smaller than those of the El and Ml transitions. In 
addition to the El-El and Ml-Ml processes, the El-Ml interference process could have 
finite contribution to the optical absorption through the mixing of the 3<i 4 4p-configuration 
to the 3(i 5 -configuration, as illustrated in Fig. [2l Such mixings are the result of the non- 
centrosymmetric environment on Fe atoms. In order to describe such processes, we employ 
a cluster model of FeC>6, which includes all the 3d and 4p orbitals of Fe atoms and the 2p 
orbitals of O atoms. The Coulomb interaction and the spin-orbit interaction are taken into 
account in the 3d orbitals. Since Fe atoms are located in the noncentrosymmetric environ- 
ment, the Ap and 3d states could be coupled to each other. A similar cluster model has been 
considered in the analysis of resonant x-ray scattering in magnetite,— where Fe atoms at A 
sites are in the noncentrosymmetric environment, at the center of tetrahedrons of O atoms. 
Deriving an effective hybridization between the 4p and 3d states as well as a ligand field on 
the 3d states through the hybridization with the O 2p states, we diagonalize the Hamiltonian 
matrix in the 3<i 5 - and 3<i 4 -configurations to obtain the energy eigenstates. These states are 
used to calculate the absorption spectra. 

In the experiment, the magnetic field was applied along the ±c axis, and the difference 
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FIG. 2: Illustration of the El-Ml interference process. The black circle indicates the presence of 
an electron in the 4p band. 

of the absorption spectra between the two directions was measured, which would be termed 
as "magnetoelectric" spectra.— Since the compound is a ferrimanget, reversing direction of 
applied magnetic field results in reversing the direction of the local magnetic moment on Fe 
atoms. Neglecting a small deviation from a perfect antiferromagnet, we simply assume that 
the direction of the local magnetic moment is simply reversed. We derive a formula for the 
magnetoelectric spectra which arise from the El-Ml process. Using this formula, we discuss 
various symmetry relations for the El-Ml process and the relation to the nonreciprocal 
directional dichroism and the anapole moment on these bases. Finally, we carry out a 
microscopic calculation of the spectra arising from the El-Ml process using the results of 
the FeC>6 cluster model. We find the spectra as a function of photon energy in agreement 
with the experiment.— 

This paper is organized as follows. In Sec. II, we introduce a cluster model around Fe 
atoms. In Sec. Ill, we describe the optical transition operators associated with Fe atoms. In 
Sec. IV, we derive the formulas of the optical absorption, and present the calculated spectra 
in comparison with the experiment. The last section is devoted to concluding remarks. 

II. ELECTRONIC STRUCTURES AROUND Fe ATOMS 
A. Crystal electric field 

We start by examining the crystal electric field around the off-center position O' = (0, 0, 5) 
to see the effect of lowering symmetry from the cubic to trigonal ones. Let charge q be placed 
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at the apexes of the octahedron. Then, the electrostatic potential 4>(x, y, z) is expanded as 

<j ) (x,y,z) = V + 5V 1 + 5 2 V 2 + ---, (2.1) 



with 

V = — - -^Xr{35z 4 - 30zV + 3r 4 ± 20V2z(x 3 - 3xy 2 )}, (2.2) 

V i = - vf{^ 3 - 3zr " ± - 3a; ^)}' ( 2 - 3 ) 

6r Q 4 

V 2 = -J{2z 2 -(* 2 + , 2 )}, (2.4) 
'o 

where the x, y and z axes are along the crystal c, a, and b axes, respectively, with the origin 
O'. The distance between the center of the octahedron and the apexes is defined as r and 
r = ^ x 2 + y 2 + z 2 . The upper and lower signs correspond to the octahedron on the left 
and right panels in Fig. HJ respectively. Term V represents the so-called cubic field term, 
which gives rise to a splitting of energy between e g and t 2g states in 3d orbitals. Term V\ 
gives rise to a coupling between 3d and 4p states, and V 2 gives rise to additional splittings of 
energy within the 3d states as well as the 4p states. These forms are inferred to be correct 
in symmetry point of view, but the covalency between Fe and O is, however, expected to 
give rise to a similar but much larger effect. We neglect the small point charge effect, and 
consider only the covalency effect discussed in the following. 



B. Hamiltonian for a FeC>6 cluster 

We now introduce the Hamiltonian of a FeC>6 cluster, and derive the ligand field on the 
3d states and the effective hybridization between the 3d and 4p states. With the 2p states 
in O atoms in addition to the 3d and 4p states in the Fe atom, we write the Hamiltonian as 

H = H 3d + H 2p + H™~ 2p + H 4p + H^~ 2p , (2.5) 



where 



H 3d = ^ E d m S ma d ma + - g {y x v% z/ 3 z/ 4 ) S vx d\^d Vz 

+ CM ( ma \ L -S\ m ' a ') d ma d rnia'-+H X c-^(S)aa>dl na d rn(7 ,, (2.6) 



mm'aa' 
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H 2 P 



(2.7) 



■3d-2p 
hyb 



^ ] ^mr] P (j)d'l n aPjr]a + H.C., 



(2.8) 



e 4p (k)pJ, ()/(T p kr; , ., 



(2.9) 



kr/<r 



-4p-2p 
hyb 



5^ ^ 2p ^KW + Hx - 



(2.10) 



The H 3 describes the energy of 3d electrons, where d ma represents an annihilation operator 
of a 3d electron with spin a and orbital m (= x 2 — y 2 ,3z 2 — r 2 ,yz, zx,xy). The second 
term in Eq. (12.61) represents the intra-atomic Coulomb interaction with the matrix element 
g [y\v 2 \ y^vi) expressed in terms of the Slater integrals F°, F 2 , and F 4 (v stands for (m, a)). 
The third term in Eq. (12.61) represents the spin-orbit interaction for 3d electrons. We evaluate 
atomic values of F 2 , F 4 , within the Hartree-Fock (HF) approximation,— and multiply 
0.8 to these atomic values in order to take account of the slight screening effect. On the other 
hand, we multiply 0.25 to the atomic value for F°, since F° is known to be considerably 
screened by solid-state effects. The last term in Eq. (12.60 describes the energy arising from 
the exchange interaction with neighboring Fe atoms, where (S)^ represents the matrix 
element of the spin operator of 3d electrons. The exchange field H xc here has a dimension 
of energy, and is ~ ksT c /4 with T c ~ 250 K. Note that this term is served as selecting the 
ground state by lifting the degeneracy and therefore the spectra depend little on its absolute 
value. The H xc is directed to the negative direction of the c axis at Fel sites when the 
external magnetic field is applied along the positive direction of the c axis. 

The H 2p represents the energy of oxygen 2p electrons, where Pj va is the annihilation 
operator of the 2p state with rj = x,y,z and spin a at the oxygen site j. The Coulomb 
interaction is neglected in oxygen 2p states. The H^ 2p denotes the hybridization energy 
between the 3d and 2p states. The energy of the 2p level relative to the 3d levels is determined 
from the charge-transfer energy A defined by A = E d - E p + 15U(3d 6 ) - 10£/(3d 5 ) with E d 
being an average of E^. Here U(3d 6 ) and U(3d 5 ) are the multiplet-averaged d-d Coulomb 
interaction in the 3d 6 and 3d 5 configurations, which are defined by U = F° — (2/63) F 2 — 



The H 4p represents the energy of the Ap states, where is the annihilation operator 
of the 4p state with momentum k, r] 1 = x, y, z, and spin a. The 4p states form an energy 



(2/63) F 4 . 
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TABLE I: Parameter values for a FeC>6 cluster in the 3d 5 configuration, in units of eV. The Slater- 
Koster two-center integrals are defined for the Fe atom at the center of the octahedron. 



F°(3d,3d) 6.39 (pda) 2pM -1.9 

F 2 (3d,3d) 9.64 (pdn) 2Pl 3d 0.82 

F 4 (3d,3(i) 6.03 {vwhpAp 3.5 

Csd 0.059 (ppvrWp -1-0 
A 3.3 



band e 4p (k). The density of states (DOS) of the 4p band is inferred from the K-edge 
absorption spectra^ as shown in Fig. El The H^ 2p represents the hybridization between 
the 4p and oxygen 2p states, where the annihilation operator of the local 4p orbital p' nl(7 may 
be expressed as = (1/y/No) J^kPk^'a (-^o is the discretized number of k-points). 

0.05 I . 1 . 1 . 1 




Energy [eV] 

FIG. 3: Density of states of the 4p band. It is constructed from the experimental .fT-edge absorption 
spectra^ with cutting off the low-energy tail coming from the life-time width of the core hole. The 
high-energy side is arbitrarily cut-off. The integrated value is normalized to unity. 

The hybridization matrices t^~ 2p (j) and tj~ 2p (j) are defined for the Fe atom at the 
off-center position. We evaluate these values by modifying the Slater-Koster two-center 
integrals for the Fe atom at the central position of the octahedron with the assumption that 
(pda)2 P ,3d, (pd7i)2 P ,3d oc d~ 4 , and {ppa)i P ^ pi {pp^)ip,2p oc d~ 2 for d being the Fe-0 distance.— 
Table H lists the parameter values used in this paper, which are consistent with the values 
in previous calculations for FesC^.— & 
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C. Ligand field and effective hybridization between 4p and 3d states 



Instead of directly treating H 3d ~ 2p and H 4p ~ 2p , we introduce the effective Hamiltonian to 
include the covalency effect. The ligand field Hamiltonian on the 3d states is given by the 
second-order perturbation as 

H 3d ~ 3d = ^J'dLdm'a + H.c, (2.11) 

mm' a 

with 

t™J d = £ C*^*0')/A, (2.12) 
jv 

where the sum over j is taken on neighboring O sites, and A = 3.3 eV is the charge transfer 
energy defined above. In addition to the ligand field corresponding to the cubic symmetry, 
we have a field proportional to 6 2 , which causes extra splittings of 3d levels in conformity 
with the form of Eq. (I2.1R . 

The effective hybridization between the 4p and 3d states is similarly given as 

H Ap ~ M = £ tX 3d 4j™ + H.c, (2.13) 

rj'mcr 

with 

K>™ d = E <sr*tf)C*tf)/(^ - E2p) ' (2 - 14) 

where E 4p is the average of the 4p-band energy, which is estimated as E Ap — E 2p ~ 17 eV. 
The coefficient t^ p ~ 3d is nearly proportional to the shift 5 of the Fe atom from the center of 
the octahedron, again in conformity with the form of Eq. (12. IB . 

III. ABSORPTION PROCESS ON Fe 

The interaction between the electromagnetic wave and electrons is described by 

H ™t = ~ /j(r)-A(r)d 3 r, (3.1) 

where j represents the current density operator, and the electromagnetic field A(r) for linear 
polarization is defined as 



A W = E\/^ ec q eiq - r + H - c -' ( 3 - 2 ) 



q 



with c q and e being the annihilation operator of photon and the unit vector of polarization, 
respectively. We approximate this expression into a sum of the contributions from each Fe 
atom: 

-^j(q,z).A(q,i) + H.c, (3.3) 



H- m t = 

c 
q 



with 



j(q,i) = E / e iq - (r - ri) j„„'(r-r,)d 3 (r 



al(i)a nl (i), (3.4) 



A(q,i) = J^ecjtf** (3.5) 
where the local current operator may be described by 

j nn ,(r - r<) = — [(V0*)0„/ - 0*V0 n /] A0*0 n / 

+ — cV x [#S0 n ,]. (3.6) 
mc 

The integration in Eq. (13 A\\ is carried out around site i, and a n {i) is the annihilation operator 
of electron with the local orbital with the wave function n (r — r$). The e and m are the 
charge and the mass of electron, and hS is the spin operator of electron. The second term 
in Eq. (13.61) . which describes the scattering of photon, will be neglected in the following 
discussion. The approximation made by taking account of the process only on Fe atoms may 
be justified at the core-level spectra, but less accurate in the optical region. The spectra 
arising from the magnetoelectric effect, however, are expected to be described rather well 
by the present approximation, since such effects mainly take place on Fe atoms. 

For later convenience, we write the interaction between the matter and the photon in a 
form, 

q V q i 

To be specific in connection with the experimental set-up,— we consider the situation that 
the photon propagates along the a-axis with linear polarization, as illustrated in Fig. 0J 



A. El transition 

The transition operator T(q, e,i) for the El transition is given by putting e iq '^ r_r ^ = 
1 in Eq. (13 AL Therefore it is independent of the propagation direction of photon. For 



y 




FIG. 4: Geometry of absorption. Light propagates along the a axis with polarization along the b 
axis or the c axis. The electric dipole moment is along the b axis. The sublattice magnetization is 
directed to the negative direction of the c axis at Fel sites and to the reverse at Fe2 sites, when 
the external magnetic field is applied to the positive direction of the c axis. When the external 
magnetic field is reversed, the sublattice magnetization is reversed. 



the polarization along the z-axis, the first term in Eq. (13.61) is rewritten by employing the 
following relation 

d 777 f 

' v cl 3 r = -^(en - e n r) / <j>* n z(t> n , d 3 r, (3.8) 



n dz Yn h?' 

where e n and e n i are energy eigenvalues with <p n and n /, respectively. The 4p and 3c? states 
are assigned to n and (f> n /, respectively. Hence the transition operator T E1 is expressed as 

T sl (q,e,*) = iB E ' J2 *Sl&(0<U(0 - <C(0pU01, (3-9) 
where % runs over Fe sites. The N E ^s are given by N E ^ X = 1/VE, N E y Z = l/y/E, N E ^ z2 _ r2 = 



El 



2/y/TE for the polarization along the z axis, N E ^ 2 _ y2 = 1/V5, N E ^ z2 _ r2 = -1/V15, N y>Xi , 
l/y/5, N E ZX = l/v5 for the polarization along the x axis, and N E xy = 1/a/5, N E x2 _ y2 = 
— 1/Vo, N E ^ z 2_ r2 = 1/yTE, N E yz = l/v5 for the polarization along the y axis, respectively. 
The coefficient B El is defined by 

poo 

B E1 = (e 4p - e 3d ) / r 3 J R 4p (r) J R 3 d(r)dr, (3.10) 
Jo 

where i?3d(r), -R 4p (r) are radial wave-functions of 3rf, 4p states with energy €3^, e 4p in the 
Fe atom. The energy difference e 4p — e^d is not directly related to the absorbed photon 
energy. Within the HF approximation in the ls 2 3d 5 4p 001 -configuration of an Fe atom,— 
we estimate it as B E1 « 7.7 x 10~ 8 cm-eV. 
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B. E2 transition 



The transition operator for the E2 transition is given from the second term in the expan- 
sion e iq '( r-ri ) fa 1 + iq ■ (r — ri) + ■ • • in Eq. ( 13.41) . Let the photon be propagating along the 
y-axis with the polarization parallel to the z-axis. Then we could derive a relation, 



<j>* n y— n ,d d r = -— (e„ - e„0 y 0;^-0 n ,d d r + - y ^L^d^r, (3.11) 

where hL x is the orbital angular momentum operator. The last term should be moved into 
the terms of the Ml transition. In the first term of Eq. (13.111) . the relevant states for n 
and (j) n > are both 3d states, and e n — e n t may be an order of the ligand field energy, which is 
less than 1 eV. Since (r 2 ) is estimated within the HF approximation as^ 

I r 4 R 2 3d (r)dr = 3.3 x l(T 17 cm 2 , (3.12) 
Jo 

we notice that the contribution from the E2 transition is smaller than that from the Ml 
transition discussed in the next subsection. 



C. Ml transition 

From the third term in Eq. f 1 3 . 6 f) . we have a relation 

J e^-^V x (0;S(/v)d 3 r = -iq x j ^S^e^^^r fa -iq x J 0;S0 n ,d 3 r. (3.13) 

Adding the contribution of the last term of Eq. (13.111) . we have a factor L + 2S in the 
transition operator. The 3d states are assigned to n and <j) n i. Hence the transition operator 
for the Ml transition is given by 

T M1 (q,e,z)=z\q\B M1 £ ^ mV dL(0^(<). (3-14) 

imm'acr' 

where B M1 = h 2 /2m = 3.8 x lCT 16 cm 2 • eV. For the photon propagating along the y axis 
with polarizations along the z and x axes, we have N^ m , a , = (ma\L x + 2S x \m'a') and 
N ma, m <a> = ( mcr \ ~ { L z + 2S z )\m! a') , respectively. 
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IV. CALCULATION OF ABSORPTION SPECTRA 



Restricting the processes only on Fe atoms, we sum up cross sections at Fe sites to obtain 
the absorption intensity I(ujq,e). Dividing it by the incident flux c/V, we have 

A 2 2 i 

JK, e) oc — E I (*/W|r(q, e, z)|* 9 (i))| 2 5(^ q + - (4.1) 

where T(q, e, z) = T £1 (q, e, z) + T An (q, e, z), and and \^ff(i)) represent the ground 

and the final states with energy E g and Ef at site i, respectively. The sum over / is taken 
over all the excited state at Fe sites. 

We first calculate the energy eigenstates |$ n (c? 5 )) with eigenenergy E n (d 5 ) in the 3d 5 - 
configuration, and |$ n (<i 4 )) with eigenenergy E n (d 4 ) in the 3<i 4 -configuration, by diagonal- 
izing the Hamiltonian H 3d + H 3d ~ 3d . As already stated in Sec. II, the exchange field H X c 
in Eq. (12.61) is assumed to be directed to the negative direction of the c(x) axis at Fel sites 
and the reverse direction at Fe2 sites when the external magnetic field is applied to the 
positive direction of the c axis. All the directions could be reversed by reversing the external 
magnetic field, since the actual compound is a ferrimagnet. The shift 5 of Fe atoms along 
the 6-axis is assumed 5 = 0.26A at Fel sites and 5 = — 0.11A at Fe2 sites, respectively. 

As regards the lowest energy state |<3> 9 (<i 5 )), we have the state & Ai under the trigonal 
crystal field, if we disregard the exchange field and the spin-orbit interaction. The inclusion 
of these interactions could induce the orbital moment (L x ), but its absolute value is given 
less than 0.004. Two types of octahedrons give the same angular momentum. 

Within the first order perturbation with the effective hybridization H ip ~ 3d , we could 
express the ground state \^ g (i)) and the optical final states as 

|*,(i)> = \%(d 5 )) 

+ Ei*-(^).^>^ ) _ ( ^ )+ ^ (k)) 

x (® n (d 4 ),k W \H 4p - sd \$ g (d 5 )), (4.2) 



!*/(<)> = |$/(rf 5 )) 



(rj4) + e 4p (k)) 



x ($ n (d 4 ),k W \H 4p - 3d \$ f (d 5 )), (4.3) 



with Eg = E g (d 5 ) and Ef = Ef(d 5 ) being the lowest and excited energies in the d 5 configura- 
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tions, respectively. Here |$ n (c? 4 ), krja) represents the state of four electrons in the 3d states 
and one electron on the 4p states specified by r](= x,y,z), spin a, and momentum k. The 
sum over k may be replaced by the integral with the 4p DOS. The explicit dependence on 
site % is abbreviated in the right hand side of Eqs. (14.21) and ( 14.31) . From these wave- functions 
we obtain the expressions of optical transition amplitudes at site i by 

M si (q,e,*;/) = (* / (i)|T B1 (q J e,<)|*p(<)> 

= ^<<Mrf 5 )|T E1 (q,e,z)|<MAkw> 

riWr\a 

+ ^($ / (d B )|^- 3d |$ n (d 4 ),k»7(7) 

M M1 (q,e,*;/) = (* f (i)\T M1 (q,e,i)\* 9 (i)) 

= ($ f (d 5 )\T M \q,e,i)\<l> g (d 5 )). (4.5) 

With these amplitudes, we have 

J(u; q , e) oc -J- ]T X] I m£1 (Q> e, i; /) + M M1 (q, e, i; /)| 2 5(^ q + £ 9 (d 5 ) - ^(d 5 )). (4.6) 
q i / 

Now we examine the symmetry relation of the amplitudes. First, let the propagating 
direction of photon be reversed with keeping other conditions. The magnetic field associated 
with the photon is reversed, N M1, s in Eq. (13.141) change their signs. Since other conditions 
are the same, we have the new amplitudes (M E1 )' = M E1 , (M M1 )' = —M M1 . Second, 
let the local magnetic moment at each Fe atom be reversed with keeping the same shifts 
from the center of octahedron. The reversing of the local magnetic moment corresponds 
to taking the complex conjugate of wave functions. Considering Eq. (14.41) together with 
Eq. (El), we have (M E1 )' = -(M El )* . Also, considering Eq. (0~5D together with Eq. pTTD . 
we have (M M1 )' = (M M1 )*. Third, let the shifts of Fe atoms from the center of octahedron 
be reversed with keeping the same local magnetic moment, which means the reversal of 
the direction of the local electric dipole moment. This operation gives rise to reversing the 
sign of H 4 P~ 3d but no change in the 3d states with the 3d 5 - and 3c? 4 -configurations, because 
the ligand field H 3d ~ 3d changes according to 5 2 . As a result, we have the new amplitude 
(M E1 )' = —M E1 from Eq. (H3D but no change (M M1 )' = M Ml . 
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As already stated, the direction of the local magnetic moment could be reversed by 
reversing the direction of the applied magnetic field, since the actual material is a ferrimagnet 
with slightly deviating from a perfect antiferromagnet. We define AJ(a> q , e) by the difference 
between the absorption intensity with the applied magnetic field along the positive direction 
of the c axis and that with the field along the reverse direction. From the second symmetry 
relation mentioned above, we have 



AJ(w„e) oc l^[{l E1 (q,e, ! |/)fM M1 (q ) e )! ;/) 



2 

7 « i f 

+ {M M1 (q, e, »; /)}* M E1 (q, e, <; /)] 5(^ q + E g (d 5 ) - E f (d 5 )). (4.7) 



Considering the sign change, we infer from the above symmetry relations that 

AI(w q , e) oc ± ■ p ioc(«) x M loc (i), (4.8) 

i 

where Pi oc («) and M loc (z) are the electric and the magnetic dipole moment of Fe atom at 
site i, respectively (Pi oc (i) oc Si = (0, 0, 5)). This relation may be regarded as a lowest order 
expansion with respect to <5, and M loc (z). The right hand side of Eq. (14.81) is the sum of the 
local toroidal moment t{i) (= 5i x Mi oc (i)).— 

Figure [5] shows the calculated A/(u; q , e) as a function of u; q , in comparison with the exper- 
iment. We have replaced the 5-function S(x) in Eq. (14. 7p by a Lorentzian form (7 /tc)/(x 2 +'J 2 ) 
with 7 = 0.1 eV. The calculated peak height at ~ 1.2 eV is set to be the same as the exper- 
imental one for the polarization e along the b axis. We have a two-peak structure around 
Hujq = 1.0 — 1.5 eV in agreement with the experiment, but could not reproduce a dip found 
experimentally around hw^ = 1.7 — 2.3 eV. On the other hand, without further adjustment, 
we have a considerable dip around ftwq = 2.0 — 2.5 eV for e along the c axis, in agreement 
with the experiments. 

Fe atoms are under the cubic symmetry without displacement, and the lowest and low- 
lying excited states are characterized as 6 A 1? 2 T 2 , 4 T l5 4 T 2 with neglecting the spin-orbit 
interaction and the exchange field.~ The excitation energies for 2 T 2 , 4 7\, and 4 T 2 are es- 
timated 1.34, 1.59, 2.45 eV, respectively, within the present cluster model. Note that the 
direct absorption processes & A\ — > 2 T 2 , 6 Ai — > 4 Ti, and & A\ — > 4 T 2 are forbidden. The 
displacement of the Fe atom generates a trigonal field and makes the energy levels of the 
excited states split. The spin-orbit interaction and the exchange field further modify these 
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Photon energy [eV] 



FIG. 5: Difference of absorption intensities AJ(u;q, e) as a function of photon energy Hujq between 
the applied magnetic field along the positive and negative directions of the c axis. Photon propagates 
along the positive a axis with polarization vector e along the b and c axes, respectively. Experimental 
data are taken from Ref. [15j |. 

states. The magnetoelectric spectra around 1.0 — 1.5 eV and around 1.7 — 2.3 eV might 
be interpreted as transitions to the states dispersed from 2 T 2 and 4 T 1? and those from 4 T 2 , 
respectively. 

V. CONCLUDING REMARKS 

We have studied the magnetoelectric effects on the optical absorption spectra in a polar 
ferrimagnet GaFeC>3. We have considered the El, E2, and Ml processes on Fe atoms, 
and have performed a microscopic calculation of the magnetoelectric spectra using a cluster 
model of FeC>6. The cluster consists of an octahedron of O atoms and an Fe atom displaced 
from the center of octahedron. We have disregarded additional small distortions of the 
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octrahedron. Due to the noncentrosymmetric environment on the Fe atom, we have an 
effective hybridization between the 4p and 3d states through the O 2p states and thereby the 
mixing of the 3d 4 4p-configuration to the 3c/ 5 -configuration. This mixing makes the El-Ml 
interference process survive and gives rise to the magnetoelectric spectra. We have evaluated 
the El-Ml process by using the energy eigenstates given in the 3(i 4 4p-configuration and the 
3<i 5 -configuration. The Coulomb interaction between 3d electrons and the hybridization are 
assumed to be nearly the same as previous cluster calculations.— 1 ^ We have obtained the 
magnetoelectric spectra as a function of photon energy in the optical region 1.0 — 2.5 eV, in 
agreement with the experiment. 

In the experiment, the conventional absorption spectra, a part independent of the direc- 
tion of magnetization, were measured with intensity about three orders of magnitude larger 
than the magnetoelectric part.— On the other hand, in the present approach considering 
only the local process on Fe atoms, the "total" intensity, which is given by the El-El and 
Ml-Ml processes, is estimated as merely one order of magnitude larger than that of the 
El-Ml process. This suggests that other processes such as the transition from the valence 
band to the conduction band involving Ga and O atoms may add larger contributions. As 
far as the magnetoelectric spectra are concerned, however, the present approach considering 
only the local process on Fe atoms is expected to work well, since the El-Ml interference 
process could take place only on Fe atoms. Finally, from a different point of view, we would 
like to comment that the approach of considering the multiple scattering of a 4p electron in 
the noncentrosymmetric potential and the Coulomb interaction between the 4pd A and the 
d 5 configurations may improve the above situation. The critical study is left in future. 

We have concentrated on the spectra in the optical region. In the x-ray region, the 
magnetoelectric effects have also been studied.— Since the core electron is excited 
there, the local approach in this paper would be better applicable to the x-ray region than 
to the optical region, where the E1-E2 (not El-Ml) interference process gives rise to the 
magnetoelectric spectra. It may be interesting to analyze microscopically the nonreciprocal 
directional dichroism observed in the Fe pre-if-edge x-ray absorption in GaFeO^ by using 
a similar cluster model. In this context, we would like to comment that the magnetoelectric 
effect on the resonant x-ray scattering spectra has been analyzed at the Fe pre-K-edge in 
Fe 3 04,— where Fe atoms at A sites are located at the center of tetrahedron in noncentrosym- 
metric environment. 
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